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Abstract 

Higher-order WKB methods are used to investigate the border between the 
solvable and insolvable portions of the spectrum of quasi-exactly solvable 
quantum-mechanical potentials. The analysis reveals scaling and factoriza- 
tion properties that are central to quasi-exact solvability. These two proper- 
ties define a new class of semiclassically quasi-exactly solvable potentials. 
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Typeset using REVT^X 



Quantum-mechanical potentials are said to be quasi- exactly solvable (QES) if it is possible 
to find a finite portion of the energy spectrum and associated eigenfunctions exactly and in 
closed form |1|]. QES potentials depend on a parameter J; for positive integer values of J 
one can find exactly the first J eigenvalues and eigenfunctions, typically of a given parity. 
QES systems have been classified using an algebraic approach in which the Hamiltonian is 
expressed in terms of the generators of a Lie algebra This approach generalizes the 

dynamical symmetry analysis of exactly solvable quantum-mechanical systems, whose entire 
spectrum may be found in closed form by algebraic means ||. 

In this paper we use higher-order semiclassical methods to examine the boundary between 
the exactly solvable part of the spectrum and the remaining energy levels in QES systems. 
We find that the large- J asymptotic behavior of the largest exactly known energy eigenvalue 
is particularly simple for QES potentials. This study leads to a natural generalization of 
QES potentials. We discover an infinite tower of potentials; the first is exactly solvable, the 
second is QES, and the rest share some semiclassical features of QES potentials but are not 
QES. 

The simplest QES potential |J is 

V(x) = x 6 - (47 - l)x 2 . (1) 

The Schrodinger equation, —ip"(x) + [V(x) — E]ip(x) = 0, has J even-parity solutions of the 
form 

i;(x) = e- x4 / A J2c k x 2k . (2) 

k=0 

The coefficients Ck satisfy the recursion relation 

4(J - fc)c fc _i + Ec k + 2{k + l)(2fc + l)c fc+1 = 0, (3) 

where < k < J — 1 and we define c_i = cj = 0. 

The simultaneous linear equations (|^) have a nontrivial solution for cq, c\, cj_i if the 
determinant of the coefficients vanishes. For each integer J this determinant is a polynomial 
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Pj{E) of degree J in the variable E. The roots of Pj{E) are all real and are the J quasi-exact 
energy eigenvalues of the potential in Eq. (JTJ) . We have computed the first 40 polynomials; 
we list the first seven of these below: 



Pi{E) 


= E, 


P2(E) 


= -8 + E\ 


Ps(E) 


= -64E + E 3 , 


P 4 (E) 


= 2880 - 240£ 2 + E\ 


P 5 (E) 


= 47104E - 640£ 3 + E 5 , 


Pe(E) 


= -5184000 + 331456E 2 - 1400E 4 + E 6 , 


P 7 (E) 


= -130940928E + 1529856E 3 - 2688E 5 + E 7 



The roots of Pj(E) occur in positive and negative pairs, and successive sets of roots 
interlace: 



roots of Pi(E) : 


0; 


roots of P 2 (E) : 


±2V2; 


roots of P 3 (E) : 


0, ±8; 


roots of P 4 (E) : 


±3.55932, ± 15.07751; 


roots of P 5 (E) : 


0, ±9.21135, ±23.56164; 


roots of P 6 (E) : 


±4.10132, ± 16.70778, ±33.22693; 


roots of P 7 (E) : 


0, ± 10.18750, ±25.56258, ±43.94052 



Let E C (J) represent the largest root and thus the critical energy that marks the upper 
edge of the quasi-exact spectrum. A numerical fit to the large-J asymptotic behavior of 
E C {J) for 1 < J < 40 using Richardson extrapolation § gives E C (J) ~ 3.0792014355 J 3/2 as 
J — > oo. We recognize the numerical constant: 

E C (J)~^-V3J 3 / 2 (J-oo). (4) 
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One may verify this result analytically by finding the minimum of the potential V(x) = 

x 6 -{AJ-l)x 2 . Minima occur at x = ±[(4J-l)/3] 1 / 4 ; at these values V f V3J 3 / 2 (J -> 

oo). Since the quasi-exact energy eigenvalues occur in ± pairs and the zero-point energy 
is negligible for large J, the asymptotic result in Eq. (Q) is confirmed J7|. However, this 
approach is not useful for all QES systems because in general the spectrum is not symmetric 
under E — > —E. 

WKB theory provides a more general derivation of Eq. (f|) and gives higher-order correc- 
tions. Furthermore, semiclassical analysis reveals features of V(x) that are generic to QES 
potentials. The leading-order WKB quantization condition is 

(n + -)?r ~ r° dx ^E-x 6 + (4J - l)x 2 (5) 
2 J-x 

for large n, where ±xo are turning points [xq is the positive zero of E — V(x)] and n is the 
quantum number of the energy level. We seek the Jth even-parity energy level E = E C (J) 
for which n = 2J. 

To evaluate the integral on the right side of Eq. ([D for large J, we begin by scaling the 
variables E c and x. The scaling E c = J 3 ^ 2 a and x = J l ^ 4 y extracts a factor of J from the 
integral and reduces Eq. (|5]) to an exact expression for the numerical constant a: 

/Ho I 
dy sja - y G + 4y 2 , (6) 
-yo 

where yo is the positive root of a — y 6 + 4y 2 = 0. 

In general, the integral in Eq. (||) is not an elementary function of a. However, for the 
special value a = the polynomial a — y 6 + Ay 2 factors into a product of a linear term 

times a perfect square: 



H*-, + V-(^)(^)'. 

Using this factorization, we express the integral in Eq. @ in terms of Beta functions: 



(7) 



r yo , i ; — 7^ 8 f , l + 2z r — 

/ dyJa-y 6 + 4y 2 = - dz — ^Vl 3 
J-yo 3 Jo Jz 



-yo 

8/3 1\ 16 „ /3 3\ 



Thus, the WKB quantization condition ([]) is satisfied with a = if \/3, which confirms the 
leading asymptotic form of the critical energy E C (J) given in Eq. (|j). 

Two aspects of the above calculation of the WKB integral are crucial. First, the QES 
potential is such that in the limit J —>■ oo, the parameter J scales out of the integrand. 
Second, the factorization in Eq. (|7|) enables us to evaluate the scaled integral and to obtain 
the factor of 7T, which then cancels from the leading-order WKB quantization condition. In 
fact, merely demanding that a — y 6 + 4y 2 factor into the product of a linear term in y 2 times 
a square of a linear term in y 2 uniquely specifies the value of a, which in turn determines 
the large- J asymptotic behavior of the critical energy E C (J). 

The factorization property extends to general QES x 6 potentials. Consider the scaled 
factored form a — V(y) = (a — y 2 )(b+cy 2 ) 2 . (The conditions a > 0, b, c > ensure that there 
is just one pair of turning points and that the WKB integrals give rise to Beta functions of 
half- integer arguments.) The leading-order WKB condition is 

2nJ ~ J dx^E - V(x) 

= J J dyja~^y~ 2 {b + cy 2 ) = 2vr + ^), 

which determines b in terms of a and c: b = - — y. This gives the potential 

V(y) = c 2 y 6 + 2(3cy A + ((3 2 - 4c)y 2 , 

where (3 = - — ^f 2 . This is exactly the large- J scaled form of the general nonsingular QES 
x 6 potential (class VI in Ref. 0]). When {3 = and c = 1, we obtain the case treated in 
Eq. ®. 

This factorization is a universal feature of quasi-exact solvability. In general, at the upper 
end of the quasi-exact spectrum, E — V(x) factors into one term, which fixes the location 
of the two WKB turning points, multiplied by a second term, which is a perfect square. 
As another example, consider the QES potential V(x) = sinh 2 x — (2 J — 1) coshx. For this 
potential there are J even-parity QES states, all of the form ijj(x) = e _cosh:r P(cosha;), where 
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P is a polynomial. To find the critical energy E c at the upper boundary of the quasi-exact 
spectrum, we apply the WKB quantization condition 



/ ° dx \JE C - sinh 2 (x) + (2 J - 1) cosh(a;) 

J-XO 

~(2J+1/2)tt (J->oo). (9) 

The substitution x — 2t reveals the factorization property of the integrand in Eq. (|9[). The 
expression E c — V(x) factors into one term 4 J — 2 — 4 cosh 2 1 (this term determines the 
turning points) multiplied by another term in the form of a perfect square, sinh 2 t. This 
factorization fixes the dependence of E c on J: E c — 1 — 2J. With this factorization, the 



right side of Eq. @ simplifies to 4 /q° dt sinh ty 4 J — 2 — 4 cosh 2 £ ~ 27r J, for large J. Once 
again, ir divides out of the leading-order WKB quantization condition. 

The remarkable connection between WKB and quasi-exact solvability persists to higher 
order in WKB theory. The factor of tt cancels from the WKB quantization condition to 
all orders leaving a purely algebraic series. To illustrate this, we do a fifth-order WKB 
calculation of E C (J) for the x G potential in Eq. ([]]). We begin by scaling the energy in 
the Schrodinger equation using E C (J) = 2[(4J — l)/3] 3//2 7. After scaling the independent 
variable, we obtain the Schrodinger equation 

e 2 tfj"(u) = (4m 6 - 3u 2 - t)V(w), (10) 

where the small parameter is e = 3/ [2 (4 J — 1)]. The series representation for the scaled 
energy 7 is 

00 

7=E«n e " ; (11) 
n=0 

where 00 = 1 as a result of the factorization in Eq. (|?|). 
The WKB quantization condition to order e 5 is 

2J + -) tt = - I duQ 1 ' 2 - — I du Q " 



2 J 2eJ ^ 96 J Q 3 / 2 
e 3 r 2Q""Q-7(Q^ 2 



du v ' + . . . , (12) 



3072 J Q 7 /2 
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where Q(u) = 7 + 3m 2 — 4-u 6 and where the contours encircle the two turning points (it is 
crucial that there be only two) and the branch cut joining them. 

Next, we expand Eq. (|T2| ) in powers of e and perform the resulting contour integrals. 
Each of these integrals can be evaluated in closed form because of the leading-order fac- 
torization f(u) = —4m 6 + 3m 2 + 1 = (1 — u 2 )(2u 2 + l) 2 . Indeed, all integrals of the form 
§ duu 2rn [f(u)] l l 2 ~ n where m and n are nonnegative integers can be expressed as simple 
multiples of 7r. For example, 

j> du [f(u)} 1/2 = y, jdu[f{u)r' 2 = -^ 
du[f{u)\-^ = -^=, lduuV(u)]^ 2 - — 



9^3' J ' l8y/3' 

We evaluate the appropriate contour integrals and then use Eq. ( |T2"D to determine the 
first five expansion coefficients a n for the scaled energy 7 in Eq. flTTl) . The expansion for the 
critical energy E C (J) is 

3^ 



4J-1 



35 5^ 



+ 



8(4J-1) 2 3(4J-1) 3 
23281 88945^ 



(13) 



3456(4J- 1) 4 5184(4J-1) 5 
Numerical results are excellent; when J = 40 the exact value of E C (J) is 746.606715392..., 
while the fifth-order WKB result in Eq. (|13D gives 746.606715384... . 

The factorizability property of QES potentials has the consequence that semiclassical 
analysis can be performed to all orders. All integrals give simple multiples of 7r and no tran- 
scendental functions ever appear. This is in stark contrast with x 2N anharmonic oscillator 
potentials for which WKB analysis leads to elliptic functions |T0| ]. 

Although factorizability is a property of QES potentials, it extends to a larger class of 
potentials that are not QES. We call such potentials semiclassically quasi- exactly solvable 
(SQES). Consider the potential 

V(x) = x w - 4Jx* - 5J 4/3 x 2 . (14) 



For fixed S this potential is not QES in the conventional group-theoretic sense one 
cannot find more than one exact algebraically determined eigenstate and corresponding 
eigenvalue. Nevertheless, for large J we can impose the scaling and factorizability conditions. 
Generalizing from the factorization of the x 6 potential in Eq. (|7|), we demand that the 
quantity E — V(x) factor into a product of a linear term in x 2 multiplied by a perfect 
square. This requirement fixes the numerical value of 5 = 3(4/5) 1 ^ 3 and the asymptotic 



behavior of the critical energy E C (J) ~ |(4J/5) 5/3 (J -> oo), and we have 

in a T a 15 f^J\^ o 
- x 10 + 4 Jx A + 1 1 - 



94J3 
2~T 



y 



4\ 3 



3 /4\ 3 



(15) 



where x = J l ^y. Because the squared factor has no real zeros, there are exactly two real 
turning points whose locations are determined by the linear factor in y 2 . 

This factorization enables us to evaluate the large-J leading-order WKB quantization 
integral exactly: 



2-kJ ~ r° dx JE - V(x) 

J-XO 



-B (- -\ +-b(-- 

5 V2'2/ 5 V2'2 



5 V2' 2 



J. 



As in the case of the QES potential in Eq. (]]]), the higher-order WKB quantization integrals 
can be done exactly. 

The potential in Eq. ( |T4D belongs to an infinite hierarchy of SQES x ik+2 potentials that 
exhibit factorization and scaling properties but which are not QES for k = 2, 3, 4, ... . The 



kth such potential is 112 

V(x 



16(fc + l) 2 
{2k + l) 2 



fc+i 



2 F i\ ^> - -x ~ k ~, -x - k; 



T(k + 3/2) 
2^T(k + 2] 

T(fc + 3/2) 



2k + l 
J k+1 



2 1 2 2 
where 2 -Pi is a hypergeometric function [O]. 



2y^T(k + 2)J 







2 


fc+i N 




1 








1 2 







s 



To 



demonstrate 



the 



factorization 



property 



we 



substitute y 2 = {T(k + 3/2)/[2y/nT(k + 2) J]} 1/(fc+1) x 2 and express E - V(x) as a linear 



factor in y 2 multiplied by a square of a polynomial of degree k in y 



2. 



16(fc+ l) 2 
{2k + l) 2 



r(fc + 3/2) 

2^T(k + 2) 



k + l 



2fc+l 



j— _ v(x) 



2J^T{k + 2) 
T(A; + 3/2) 



2k + l 
fc+1 



, 2 



.n=0 



Ti!\/7r 



-?7 



Note that there are two real turning points at y — ±1. 

This factorization fixes the large-J asymptotic behavior of the critical energy: 



E c 



16(fc + 1) 



T(k + 3/2) 
2^Y(k + 2) 



fc+i 



(2k + \y 

Leading-order WKB verifies this asymptotic behavior: 
2nJ ~ r° dx - V(x) 

J-Xn 



2k + l 
J + ! . 



(16) 



2^(k + l)\ [i VT^ *T(n+ 1/2) fc _ n 
T(A; + 3/2) 7o ^ ^ 



r(* + 3/2) £j 



Let us examine the hierarchy of SQES potentials for integer [13] values of k using the 



general formula [11 



2 Fi(l/2, —k - 1/2; -k + 1/2; z) 

r(i/2 - fc) 



r(i/2) 

For k = we have £? c ~ 4J and 2 -Pi (§, — §; §; 



1/2; z) 




(di) 


" (1 _ 2 )fc+l/2- 







(17) 



\/l — z, and we obtain the exactly 

W V3JV 2 and 



solvable harmonic oscillator potential, V(x) = x . For k = 1 we have £? c • - 
2 F X f |, — |; — |; = (1 + 2z)\/l — z, and we obtain the large-J asymptotic approximation 
to the x 6 QES potential in Eq. (@): V(x) = x 6 - 4Jx 2 . For fc = 2 we have E c ~ §(4J/5) 5 / 3 
and 2-P1 f§, — |; — |; ^) = (l + f 2 + f-2 2 ) \/l — -2, and we obtain the x 10 SQES potential in 
Eq. (|15|) : V(a:) = x 10 — 4Jx 4 — ^(4J/5) 4 ^ 3 x 2 . In the limit fc-*oowe obtain the square- well 
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potential, for which E c ~ 4J 2 . Thus, at the bottom of the SQES hierarchy is the exactly 
solvable harmonic oscillator, followed by the QES x 6 potential, which is then followed by the 
range of new SQES potentials and, finally, in the k — > oo limit, by the square-well potential. 

For QES potentials the critical energy E C (J) lies at the upper boundary of the quasi- 
exact spectrum. For SQES potentials, the critical energy E C (J) is still a significant point in 
the energy spectrum; at this point we observe a sharp change (a first-order phase transition 
for large J) in the density of states. If we use ordinary WKB analysis to find the 2nth energy 
eigenvalue for a x 4k+2 potential, where n is not correlated with J, we find that as n — > oo 



Observe that the numerical coefficient in this asymptotic relation is a different function of k 
from that in Eq. flT6|). For example, for k = 1 we have E c ~ 3.0792 J 3 / 2 but E 2n ~ 6.4066n 3/2 . 
Also, for k = 2 we have E c ~ 3.1024J 5 / 3 while E 2n ~ 7.4235n 5 / 3 . As k — » oo we have 
E c ~ 4J 2 while E 2n ~ 7T 2 r7, 2 . However, for the special exactly solvable case k — 0, where the 
level spacing is constant, they are the same: E c ~ 4 J and i^n ~ 4n. 

Finally, we remark that SQES models are a generalization of the Ising limit in quantum 
mechanics or quantum field theory. The Ising limit is an asymptotic balance between a mass 
term and a self-interaction term: For the potential gx 4 — m 2 x 2 we take g ~ m 2 with g large 
and get a double-well; g scales out, and we obtain the Ising limit. For SQES models there 
is an asymptotic balance for large J among all the coefficients of the potential. 
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